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Heterogeneity in the degree distribution is known to suppress global synchronization in complex 
networks of symmetrically coupled oscillators. Scale- free networks display a great deal of heterogene- 
ity, containing a few nodes, termed hubs, that are highly connected, while most nodes receive only 
a few connections. Here, we show that a group of synchronized nodes may appear in scale-free net- 
works: hubs undergo a transition to synchronization while the other nodes remain unsynchronized. 
This general phenomenon can occur even in the absence of global synchronization. Our results 
suggest that scale-free networks may have evolved to complement various levels of synchronization. 
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The last decade has witnessed a tremendous growth 
of interest in various kinds of collective dynamics in net- 
works with complex structures, ranging from physical, 
biological to social and engineering systems [TUTU] . Real- 
world complex systems have been modeled as networks 
of interacting nodes. Synchronized activities have a ma- 
jor impact on the network with important fitness con- 
sequences to all nodes and network functioning. The 
network structure exerts dramatic influence on its syn- 
chronization properties [SHS]. 

Recent studies reveal that disparate real-world net- 
worked systems share important structural features such 
as the scale-free property [TTJ [T3] . Scale-free networks 
are characterized by a high level of heterogeneity in the 
node's degree - the number of connections of a node. 
Such networks contain a few high-degree nodes, termed 
hubs, while most nodes receive only a few connections. 
The hubs serve specific purposes within their networks, 
such as regulating the information flow and providing 
resilience during attacks. They severely affect the dy- 
namical processes taking place over scale-free networks, 
particularly the emergence of global synchronized motion 

[3H5]. 

Heterogeneity in the degree distribution may lead to 
a hierarchical transition towards global synchronization, 
with hubs synchronizing first, followed by the low-degree 
nodes [Bj. In large scale-free networks, however, the 
heterogeneity inhibits global synchronization [5 . This 
turns out to be a desirable property, since in most real- 
world networks where synchronization is relevant, global 
synchronization can be related to pathological activities, 
such as epileptic seizures [7] and Parkinson's disease [5] in 
neural networks. The study of collective behavior apart 
from global synchronization is thus of substantial inter- 
est. 

In this letter, we show a general cluster synchronization 
in scale-free networks - only the hubs undergo a tran- 
sition to synchronization even in the absence of global 
synchronization. Interestingly, the very heterogeneity 
that may prevent global synchronization is the primary 
ingredient of hub synchronization. We provide condi- 
tions for the onset of hub synchronization and determine 



the persistence under small perturbations. One direct 
consequence of our theoretical analysis is that hub syn- 
chronization is both dynamically and structurally stable, 
thus, allowing the network to function in a flexible and 
robust way. 

Our approach is to introduce nonlinear dynamics on 
each node and then perform stability analysis to deter- 
mine when the hubs synchronize. From the point the 
view of stability, reasonable arguments show that the 
network dynamics acts as a small noise-like coupling. 
Hence, the linear stability of the synchronized hubs is 
maintained. Later on, in the large size limit, we provide 
a rigorous treatment on the linear stability problem. Our 
analysis is based on the new results of the theory differ- 
ential equations and spectral graph theory. 

We consider a network compose of n nodes, and label 
the nodes according to their degrees k\ < &2 < • • • < k n , 
where k\ and k n denote the minimal and maximal node 
degree, respectively. Hence, the ith. node has degree fcj. A 
scale-free network is characterized by the degree distribu- 
tion P(k), the probability that a randomly chosen node 
within the network has degree k, that follows a power-law 
P(fc) = ckr 1 , for k\ < ki < k n , where c is the normaliza- 
tion factor. The degree distribution is normalizable for 
7 > 1, and for large k n we have c sw (7 — l)kj~ . The 
mean degree (k) attains a finite limit for large k n pro- 
vided 7 > 2. We consider only connected networks with 
well defined mean degree, that is, 7 > 2. 

The dynamics of a general network of n identically 
coupled elements is described by 
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A l0 [E{ Xj )-E{x t )], (1) 



here Xi G R m is the m-dimensional vector describing the 
state of the ith node (node with degree ki), F : R m — > R m 
governs the dynamics of the individual oscillator and 
is assumed to be smooth, E : R m — > R m is the cou- 
pling function (without loss of generality assumed to be 
a constant matrix), a is the normalized overall coupling 
strength [13], and A is the adjacency matrix. A encodes 
the topological information of the network, defined as 



1 if nodes i and j are connected and Ai 







otherwise. Note that A is symmetric, and by definition 

k% — Zj j Aij ■ 

We wish to show that a group of oscillators having 
nearly the same number of connections as the main hub 
may display a synchronized motion. Consider £j = 
x n — Xi, thus, synchronization is possible between the 
nodes i and n if £j — > 0. Stability of this synchronized 
state is determined by analyzing the variational equa- 
tions governing the perturbations, which read 



& = Ki(t;a)£i +arji 



(2) 



where the matrix i^i(t; a) = [-D-F(:r n (£)) — afiiE] depends 
continuously on i, DF stands for the Jacobian matrix of 
F, fii — ki/k„ is the normalized degree, and 



Vi 



z 1 ^.. 



A»i)E((i) 



is the coupling term. 

Neglecting the coupling term rji the equations govern- 
ing the evolution of the perturbations f; and are decou- 
pled from the other perturbations and read 



& = [DF(x n (t)) + amEfo 



(3) 



We now assume that Eq. (pi) is Lyapunov regular and 
that its fundamental matrix is integrally separated |14j . 
The stability of the zero solution of Eq. ^ is determined 
by its largest Lyapunov exponent A(a//j), which can be 
regarded as the master stability function of the system 
[31 E]. The perturbation & is damped out if A(a fii) < 0. 

For many widely studied oscillatory systems the mas- 
ter stability function A(a/ij) is negative in an interval 
a\ < afii < a2 for general coupling function E [31 0]. 
The perturbation & is damped out if ct\ < a/J,i < a?. 
Moreover, normalization imposes (x n = 1 and \i\ ex fc~ , 
hence, as A; n increases, /Ui converges to zero. Not only 
/Xi, but most of the normalized degrees /ij will converge 
to zero. Therefore, it will be impossible, for large k n , to 
have o.\ < ctfii < a 2 for all i = 1, 2, • • • , n. Hence, in the 
thermodynamic limit no stable global synchronization is 
possible in scale- free networks. 

Now take a in the stability region. Then, the state 
x n = x n —\ is linearly stable. This is true as long as we 
can neglect the coupling term rji. Under the effect of 
r/i local mean field arguments show that x n ~ x n -i is 
stable. The argument goes as follows. If A(a/J, n -i) < 0, 
we guarantee the linear stability of £ n _i. Moreover, if the 
remaining oscillators are not synchronized, the coupling 
term r] n -i can be viewed as a small coupling noise, as 
long as the signals Xi are uncorrelated, with a fixed and 
k n large |16) . Results from ordinary differential equations 
state that the linear stability is maintained under small 
perturbations [121 HZ]- Therefore, if at t = we have 
x n (0) — x„_i(0) ~ 0, then for all t > it yields x n (t) — 
x n -i{t) « 0. 



These arguments cannot be applied to low-degree 
nodes. The reason is that to set the low-degree nodes 
into the stability region we must have a/xi w ai, requir- 
ing a to be as large as k n . Hence, the coupling term ar\i 
cannot be made small for low degree nodes. 

The mean field arguments also hold for correlated scale- 
free networks. The node correlation does not play a ma- 
jor role to the onset of hub synchronization. For instance, 
the Barabasi-Albert (BA) scale-free model is known to 
present finite size node correlation, hubs are likely con- 
nected [12) . If we rewire the connections between the 
hubs, connecting the hubs with the low degree nodes, 
the mean field argument is still valid, that is, hub syn- 
chronization still takes place. 

We illustrate this phenomenon with numerical exper- 
iments. We generate a Barabasi-Albert (BA) scale-free 



network with 3 x 10 nodes and 



3 12 . The network 



has largest degrees k n = fc„_i = 165. Each node x% is 
modeled as a Rossler oscillator, for x% — (xu,X2i,X3i) T 
we have F(xi) = [x 2i - x 3i ,x u + 0.2x 2 i,0.2 + x 3i (xu - 
7)) T . We consider E to be a projector in the first compo- 
nent, i.e., E(x,y,z) T = (a;,0,0) T . The master stability 
function A(a) has a stability region for a € («i, a^) with 
a\ ~ 0.13 and ct2 ~ 4.55. Global synchronization in this 
network is impossible 18J. 

For a — 0.30 we have observed the hub synchronization 
x n ~ x n -\. In Fig. [TJa) the time series x n is depicted 
in full line while x n -\ is depicted in light gray line and 
?7„_i in bold line. Fig. [I] (a) shows that the local mean 
field approximation on ?7„_i indeed holds, as shown in 
the times series x n _i ~ x n . In Fig. [TVb), the differences 
£„_i = x n -\ — X n is depicted in bold line while £2000 = 
2^2000 — %n in full line. Clearly £„_i « whereas £2000 
presents large fluctuations. 

All this reasoning can be set into a rigorous frame in 
the thermodynamic limit, for uncorrelated scale-free net- 
works. To tackle the problem let us introduce Q(t) — 
Xi(t) — s(t), where s(t) is a given typical trajectory of 
x = F(x). Consider £ = (Ci, £2, . . . , Cn) T and /1 = 
diag(ni,^2) • • • jMn)- Hence, C € M. mn . The variational 
equations of the perturbations £ can be written in a con- 
venient block form 



C = n(t:a)C + aBC 



(4) 



where tt(t; a) = I n ® DF(s(t)) — afi ® E, with ® stand- 
ing for the Kronecker product, and B = k~ l A®E is 
the coupling among the variational equations. We shall 
demonstrate that for large scale- free network with 7 > 2, 
the term coupling term can be made arbitrarily small. 

According to the aforementioned arguments ft(t; a) 
splits into independent blocks as in Eq. (pi). By choos- 
ing a fixed a such that nodes with degree larger than 
fc„_£ have their perturbations damped out, we guaran- 
tee that £ nodes display a synchronous behavior with the 
main hub x n . In other words, M. nm = U © S, where U 
and S respectively the unstable and stable spaces, clearly 
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FIG. 1: [Color online] Hub synchronization in a BA scale- 
free network of 3000 coupled Rossler oscillators with coupling 
parameter a = 0.3. a) Time series of the largest hub x n 
(full line) and the second largest x n -i (light gray line). The 
coupling term rj n -\ (bold line) spoiling the stability of the hub 
synchronization is small as predicted by the local mean field 
arguments, b) Time series of the largest hub x n (full line) and 
of a low-degree node X2000 (light gray line) The corresponding 
node degrees are k„ — k n -i = 165 and feooo — 3. 



dim(E7) = (n — t)m and dim(S') = tra. Notice that on 
the subspace S all Lyapunov exponents are negative. 

It remains to show that the coupling term can be made 
as small as one wishes whenever k n is large enough. Thus, 
results of qualitative theory of ordinary differential equa- 
tions guarantee that the linear stability is not affected by 
small continuous perturbations [19 . 

By our hypothesis on the symmetry of the matrix A 
the spectral theorem guarantees that 

A = NJN' 1 . 

where N is an orthogonal matrix and J = 
diag(Xi, A2, . . . , A n ) is the matrix of the eigenvalues of A 
ordered according to their magnitudes Ai < A2 < • • • < 
A n . 

We endow the vector space W mn with the norm || • ||* 
such that for u € M. mn we have \\u\\* = \\N <E) I m u\\oo, 
where ||u||oo = su Pi \ u i\ for i = 1,2, ••■ ,nm. We also 
make use of the induced matrix norms. Now we claim 
that given 5 > there exists K such that for all k n > K 
we have 

ll-BH* < S - 

Indeed, by using the induced matrix norm we can obtain 
bounds in terms of the largest eigenvalue of A. We post- 



pone the technical details and go directly to the result 
which reads \\A(g> E\\* < X n \\E\\oo- 

Under mild conditions [5D] the largest eigenvalue of a 
scale-free network scales almost surely as A„ = k@, where 
depends on 7. We have two distinct cases: (i) (3 = 3 — 7 
for 2 < 7 < 2.5; and (ii) j3 = 1/2 for 7 > 2.5. Putting 
all estimates together yields 



\B\ 



1 



k 1 



•P' 



(5) 



Hence, for k n large enough our claim follows. 

This analysis is grounded on the fact that X n /k n — > 0. 
This is also the case for correlated scale-free networks [10] , 
whenever the correlations preserve the scale-free charac- 
ter. These moderate correlations are immaterial for hub 
synchronization, as finite size correlation in the BA scale- 
free model. 

In summary, we analyzed a general phenomenon in the 
synchronization of large scale-free networks, namely, the 
synchronization of hubs even when the entire network 
is out of synchrony. Our theoretical analysis provides in- 
sights into further generalizations for the master stability 
function. The stability analysis of the synchronous hubs 
can be tailored to the master stability function and the 
coupling term due to the underlying network dynamics. 
We have shown that for large scale- free networks the cou- 
pling term can be controlled, effectively acting as a small 
noise-like perturbation on the hubs. 

Hub synchronization has counterintuitive effects. For 
example, the hubs do not need to be directly connected 
to synchronize. Remarkably, when the hubs synchronize, 
the low-degree nodes are out of synchrony; these nodes, 
however, are responsible for mediating the exchange of 
information between the hubs. This seems to challenge 
our understanding of the role of synchronization in the 
exchange of information within complex networks [9] . 

We believe that our findings provide strong evidence 
that incomplete, hub-driven, synchronization may be at 
least as important and persistent in real-world networks 
as other forms of synchronization and collective behaviors 
previously examined in the literature. 
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